dir = here()
dir
library(data.table)
library(ggplot2)
library(gridExtra)
library(reshape2)
library(matrixStats)
library(reshape2)
library(here)
text.size = 10
dir = here()
dir
wd = gsub("Figures", "BLS_prices", dir)
dir = here()
wd = gsub("Figures", "BLS_prices", dir)
setwd(wd)
d = fread("prices.csv")
d$name = factor(d$name)
d$name =  reorder(d$name, -d$index, mean)
p.log = ggplot() +
geom_line(data = d, aes(x  = year, y = index, col = name)) +
scale_x_continuous(breaks = seq(1900, 2060, 20)) +
scale_y_log10("Price Index (1935 = 1)") +
ggtitle("log scale") +
theme_bw() +
theme(panel.border = element_blank(),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position =  "none",
legend.margin=margin(0,0,0,-5),
plot.title = element_text(face = "bold", hjust = 0.5, vjust = -2, size = rel(0.9)),
axis.line = element_line(color = "black"),
axis.title = element_blank(),
axis.text.x = element_text(margin=margin(5,5,3,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size*0.8, family="Times"))
p = ggplot() +
geom_line(data = d, aes(x  = year, y = index, col = name)) +
scale_x_continuous("Year", breaks = seq(1900, 2050, 20)) +
scale_y_continuous("Price Index (1935 = 1)") +
ggtitle("A.  Divergent Price Change") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.text = element_text(size = rel(0.7)),
legend.key.height = unit(0.3, "cm"),
legend.key.width =  unit(0.3, "cm"),
legend.title = element_blank(),
legend.position = c(0.25, 0.6),
plot.title = element_text(face = "bold", hjust = 0.5, size = rel(0.9), vjust = -0.5),
axis.line = element_line(color = "black"),
axis.title.x= element_text(vjust=-0.4, size=rel(0.9)),
axis.title.y = element_text(vjust= 1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,3,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
wd = wd = gsub("Figures", "Vintage_gdp", dir)
g = fread("vintage_gdp.csv")
wd
setwd(wd)
g = fread("vintage_gdp.csv")
last = subset(g, g$year == max(g$year,na.rm = T) )
spread = (max(last$gdp) - min(last$gdp)) / mean(last$gdp)
uncertainty = paste("Uncertainty ==", round(spread*100), "*\'%\'")
gdp = ggplot() +
geom_path(data = g, aes(x = year, y = gdp, col = series_year), size = 0.05, alpha = 0.7) +
scale_color_gradientn("Base Year*", colors = rainbow(8)) +
scale_x_continuous("Year") +
scale_y_continuous("Real GDP (1947 = 1)", seq(1, 10, 1)) +
ggtitle("B.  Divergent Measures of US Real GDP") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.text = element_text(size = rel(0.65)),
legend.title = element_text(size = rel(0.7)),
legend.key.height = unit(0.5, "cm"),
legend.key.width =  unit(0.3, "cm"),
legend.position = c(0.2, 0.6),
plot.title = element_text(face = "bold", hjust = 0.5, size = rel(0.9), vjust = -0.5),
axis.line = element_line(color = "black"),
axis.title.x= element_text(vjust = -0.4, size=rel(0.9)),
axis.title.y = element_text(vjust= 1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,3,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times")) +
annotate("text", label = "Real GDP Growth\n1947-2018",  x = 2010, y = 4.5, family = "Times", size = 2)  +
annotate("text", label = uncertainty, parse = T, x = 2010, y = 3.7, family = "Times", size = 2)
setwd(dir)
gA =  ggplotGrob(p)
gB =  ggplotGrob(gdp)
setwd(dir)
png("prices.png", width = 3.2, height = 5.5, units = 'in', res = 600)
grid.arrange(rbind(gA, gB, size = "last"))
dev.off()
source('~/Desktop/Supplementary Material/Figures/prices.R')
library(data.table)
library(ggplot2)
library(gridExtra)
library(reshape2)
library(matrixStats)
library(reshape2)
library(here)
text.size = 10
dir = here()
# selected CPI indexes
##################################################################
wd = gsub("Figures", "Data/BLS_prices", dir)
setwd(wd)
d = fread("prices_selected.csv")
View(d)
10/0.38
library(data.table)
library(here)
dir = here()
setwd(dir)
url = fread("urls.csv")
View(url)
